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Abstract. We introduce a formalism for constructing 
models of the gravitational lens action of arbitrary ellip- 
tical surface mass distributions including varying elliptic- 
ity and isodensity-twist. The projected elliptical surface 
mass distributions are cut to slices of constant surface 
mass density for which the deflection is easy to compute. 
The deflection of the total mass distribution is found by 
superposition of the action of all slices. 
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1. Introduction 

There is an ongoing discussion on how to model elliptical 
mass distributions like elliptical galaxies and clusters as 
gravitational lenses. For these models the 2-dimensional 
potential of the projected mass distribution is needed. For 
some elliptical 3-dimensional mass distributions (as, e. g., 
triaxial, homoeoidal) the isodensities of the projected sur- 
face mass distribution are exact elliptical curves. For oth- 
ers (e. g. confocal) this is only approximately true (see 
e. g. Schneider et al. 1992 or Schramm 1990, hereafter Pa- 
per I). 

Some authors prefer to start with 2-dimensional ellip- 
tical potentials which are easy to compute (Blandford & 
Kochanek 1987, Kochanek & Blandford 1987). Their dis- 
advantage is that the related surface mass distributions 
do not share their symmetries and that they are not ap- 
plicable for mass distributions with high ellipticity. 

Other authors prefer to start with elliptical mass distri- 
butions. Bourassa et al. (1975, hereafter BK) introduced a 
complex formalism describing the lens action of projected 
spheroidal (and homoeoidal) mass distributions with ar- 
bitrary ellipticity. Schramm (Paper I) found that their re- 
sults are also valid for projected triaxial distributions and 
introduced an equivalent real formalism. Closed solutions 
have rarely been found and are by construction restricted 
to the homoeoidal symmetry, where the isodensities are 



similar and similarly situated, concentric ellipses (see Kas- 
siola & Kovner 1993, hereafter KK, for recent results or 
Narasimha et al. 1982 and below). 

Both approaches differ only slightly concerning im- 
age configurations produced by elliptical lenses of mod- 
erate ellipticity (see also KK). However, using potentials 
of elliptical surface mass distributions allows more free- 
dom in modelling elliptical galaxies and clusters. On the 
one hand some special elliptical surface mass distributions 
become tractable for theoretical purposes, on the other 
hand the isophotes of observations of elliptical galaxies 
are typically fitted by ellipses using standard packages like 
STSDAS. The results are three functions of the elliptical 
radius (introduced below), namely intensity, ellipticity and 
isophote-twist which may be used as input for modelling 
gravitational lenses (see e.g. Bernstein et al. 1993). Even 
if the mass-to- light ratio is not known, the mass might 
share the symmetry of the isophotes (compare e.g. Peletier 
1989) and therefore only the mass-to-light ratio as a func- 
tion of radius must be fitted. A failure of this fit would 
indicate different symmetries for light and mass. 

Our program is to take the 2-dimensional, projected 
surface density distribution with elliptical symmetry and 
cut it into slices of constant density (see Fig. 1). This 
means that we substitute the density distribution by a 
step function. For approximative purposes this function 
could be chosen as the lower, upper or symmetric case of 
the step function where the position of the steps are cho- 
sen as given by observations or computational needs. (In 
Fig. 1 we show the lower step function as an example.) 
For analytical purposes the choice of step function does 
not matter as long the density is Riemann integrable. We 
therefore define an elliptical slice to be a bounded ellipti- 
cal area of constant surface mass density. (We avoid the 
expression elliptical disk which could be misleading.) For 
these elliptical slices we present solutions for the lens ac- 
tion including potentials (for the time delay) and deriva- 
tives (to find the critical and caustic structure of the lens 
mapping). It is then shown how to superpose elliptical 
slices to obtain a description for arbitrary elliptical sur- 
face density distributions. If the slices are taken to be of 
finite thickness, we obtain an approximate, numerical de- 



Fig. 1. The surface mass density <t(Ye) as a function of the 
elliptical radius defined in Eq.(6) in arbitrary units. Also given 
is a related lower step function s(Ye). The shaded region rep- 
resents an elliptical slice of thickness Act 



scription. An appropriate, exact analytical description is 
found by taking the limit of infinitely thin slices and inte- 
grating over all their contributions. 

2. Basics 

The gravitational lens action of a deflector at angular size 
distance D d on a source at distance D ds from the deflector 
and D s from the observer is described by 



w(z, z) = z — V$(z, z), 



(1) 



where the complex quantities w = u + iv,z = x + iy de- 
note angular positions in the source and deflector plane, 
respectively. The bar denotes complex conjugation and the 
gradient is defined by V = d/dx + id/dy = 2d/dz. 

Although we introduced in Paper I a real formalism 
for the purpose of making clear the geometrical meaning 
of the occurring quantities, we use here the complex nota- 
tion for several reasons: (a) it is typically not necessary to 
split the lens equation into real and imaginary parts un- 
less we need numerical data, which is normally automati- 
cally achieved using complex variables in e.g. FORTRAN; 
(b) some results are found in a very simple and elegant 
form, e.g. the lens action of an elliptical slice or the para- 
metric equations for the critical curves as introduced by 
Witt (1990); (c) if it is unavoidable to look at real and 
imaginary parts for analytical purposes, the equations can 
be transformed using powerful symbolical algebra systems 



ware on an AMIGA 3000~'~'^' computer system for evalu- 
ating and checking our formulas and for preprocessing the 
figures. 

The real potential $ obeys the Poisson equation 



A$(z,z) = 2<r(z,z), 



(2) 



where the dimensionless real surface density is defined as 



<r(z,z) 



£(z,z) 



(3) 



and the critical surface mass density is as usual defined by 



c 2 D s 
4ttG D d D ds 



(4) 



where all distances Dij are measured as apparent size dis- 
tances. 

An ellipse is defined by 



b 2 



(1 + 2 



: 1 , or 

y 2 
(i-0 2 



(5) 



(6) 



respectively. The parameters a, b are the major and minor 
axis of the ellipse, which is evidently also fully determined 
by the elliptical parameter e and the elliptical radius te 
(note that the e-parameter is neither equal to the ellip- 
ticity which is 1 — b/a, nor equal to the numerical eccen- 
tricity which is \/ a 2 — b 2 /a, but is related to both). For 
convenience we give some useful relations between both 
descriptions Eqs. (5) and (6): 

(7) 

= — , (8) 
4«i. (9) 



a + b 
a + b 



and for the foci ±/ of the ellipse 



Note that the relation between / and r is restricted by 
< e < 1 so that for a constant / a lower limit r m ; n = //2 
exists. 

A surface mass distribution a = <r{r-^) with e = const, 
is called homoeoidal, and with / = const, it is called con- 
focal. The isodensities due to a confocal distribution are 
given by a set of confocal ellipses, parameterized by A, 
confocal with respect to the ellipse given by Eq.(5) 



y 
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(10) 



For given a,b we find the confocal ellipse at x,y solving 
Eq. (10) for A. 



Elliptical slices acting as gravitational lenses are interest- 
ing for several reasons. As we shall see they represent in- 
side a pure a — j (convergence-shear) lens with two critical 
densities (for a constant density: two perpendicular focal 
lines at different distances). Taking one slice alone (e.g. 
as a cluster model), it could therefore produce very elon- 
gated images of a source. A detailed discussion of this lens 
type will be published in a subsequent paper (Schramm 
& Kayser 1993). For our purposes elliptical slices will be 
the basic element of our construction set. 

We derive the lensing properties of an elliptical slice as 
a trivial consequence of the results obtained in BK. They 
could also be found following the track of the Appendix 
A in Paper I (see below). Blandford et al. (1991) obtained 
later an equivalent result for the lens action of an inclined 
disk using elliptical coordinates (compare also Appendix 
B of Paper I). 

3.1. A short summary of BK's results 

BK introduced a formula describing the lens action of a 
homoeoidal mass distribution 



I(z,z) = 2(1 - e 2 )sign(z) 



rE(z) a(r)r<\r 



4er 2 



(11) 



where r^z) is given by Eq. (6). We note that we integrate 
over all matter inside the ellipse hit by z. The outer parts 
do not contribute to the deflection. The integrand could be 
interpreted as the lens action due to an infinitesimal thin 
homoeoidal ring of constant surface mass density. For de- 
tails see the Appendix B of Paper I. A generalized signum 
function sign(z) was introduced by Bray (1984) to assure 
the right symmetry behaviour of the BK formula. The def- 
initions given below assure the equivalence to the results 
obtained with real formalism of Paper I, if the canonical 
choice for the principal value (as taken by FORTRAN77 
compilers) is chosen for the evaluation of complex square 
roots. Other choices are possible, however. Bray's defini- 
tion also ensures that the internal and external solutions 
obtained below (Eq. (15), Eq. (17)) coincide at the bound- 
ary of the slice. 

Note that the following definition differs from the 
canonical signum function as typically implemented in 
computer (algebra) languages. Bray defined: 



sign(z) 



(12) 



+ 1 if x > 
— 1 otherwise 



(x = and y ^ 0) 



Note that sign(z) changes only the sign of the t/-axis, 
that 



V$ = I 



3.2. The deflection due to an elliptical slice 



(14) 



For a = (To = const. Eq.(ll) is easily evaluated. We de- 
scribe the lens action of an elliptical slice with elliptical 
radius te- There is an external solution for z lying outside 
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(15) 
(16) 



and an internal solution for z lying inside. We note that 
we have only to integrate up to the elliptical radius corre- 
sponding to z. With z = x — iy and Eq. (6) it follows that 
(compare Eq. (A3) of KK) 



/in = ((l-e)z + i(l + e)s/)<7o 
= (z - ez) (T . 



(17) 



Note that 7; n does not depend on the actual size of the slice 
but on the elliptical parameter e defined by the bounding 
ellipse. Equations (15-17) are equivalent to the results ob- 
tained in Appendix A in Paper I. The deflection angle a 
due to an elliptical slice reads with our normalization: 



a = a x + ia v 



ab fx . y ' 



(18) 



where the primes denote the confocal semi-axes defined 
by Eq.(10). For a proof of the equivalence one has to show 
that 



a = —I 



(19) 



For the internal solution this is a trivial task since in 
Eq. (18) the primes are dropped. For the external solu- 
tion one can show after some algebra that Eq. (19) leads 
to Eq. (10). We repeat the result from Paper I since in the 
real version it is easier to see that the norm of the deflec- 
tion has homoeoidal symmetry inside and confocal sym- 
metry outside the deflecting slice. In addition it is easy to 
see that with increasing |z| the deflection angle ( in this 
confocal region) converges monotonically to the deflection 
angle of a point mass. 

3.3. Derivatives and critical structures 



sign(z)sign(z) 



-1 if 







+ 1 otherwise 



(13) 



Calculus of non-analytic complex functions could be ex- 
pressed in terms of Wirtinger-derivatives (see e.g. Janich, 
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(22) 
(23) 



For the elliptical slice these derivatives become very 
simple since dl ex /dz = outside and <97; n /<9z = <7o inside 
(compare also Witt 1990 and Schramm et al. 1993 for de- 
tails). The derivatives of the components of the deflection 
angle with respect to the coordinates (x,y) can be found 
by using Eq. (22) and splitting it in its real and imaginary 
parts. 

The determinant of the Jacobian of gravitational lens 
mappings w is given by 



dw 
dz 



dw dw 
dz dz 



(24) 



since dw/dz is always a real function. For the lens map- 
ping due to an elliptical slice we find for J in the inner 
region simply 



Jin = (1 " to) 2 - (ea ) 2 



(25) 



which shows that we find two focal lines, J = 0, as typical 
for pure astigmatic lenses, corresponding to two critical 
values for the density 



1 



Cl,2 



lie 



(26) 



For the external case we note that 7 ex in Eq. (15) is an 
(anti)analytical function solely of z (since there is no mat- 
ter outside the slice). Then the parametric form of the 
critical curve J ex = (Witt 1990) becomes 



57 

dz 



,0 < <j> < 2tt 



(27) 



which leads to a parametric equation of the external crit- 
ical curves 



e i*( e i* + i^<7 ) 



(28) 



Note that we have to take only those parts of these curves 
which lie inside (for Eq. (25)) or outside (for Eq. (28)) of 
the considered elliptical slice, respectively. 



In matter free regions the varying part of the potential is 
given by (compare also Appendix B in Paper I) 



$ ex = Re J 7 ex (z)dz 
which yields 



(29) 



-sign(z)z^ z 2 - f 1 + z 2 ^j (T 



(30) 



The constant of integration was chosen to achieve a varia- 
tion as In |z| far from the slice. Inside the slice the potential 
is given by 



Re (7 m ) dx + / Im (7 m ) dy 



which is simply 

l - ((l-e)x 2 + (l + e)y 2 +) a + C 
where the constant C is chosen as 
C7 = ,7orl(l- e 2 )ln(rE) . 



(31) 



(32) 



(33) 



With this choice, in order to give a continuous variation 
across the boundary of the slice, the potential is smooth 
and the boundary of the slice turns out to be only a point 
of inflection. As shown in Fig. (2) the potential shows inte- 
rior a homoeoidal symmetry but with lower eccentricity as 
the slice. Far from the slice it tends towards the potential 
of a point mass. 

3.5. Rotated and shifted slices 

We can express the lens action of an elliptical slice which 
is rotated by an angle <p with respect to the (x, t/)-frame. 
The relation for the lens action between the original and 
the rotated frame (x',y') reads: 



7 = e lv (I)' 



(34) 



note that in general (7)' ^ (7'). With z' = ze llp and 
(z') = ze llp we find in the external region 



7 



l-e 2 
2e 



ze 2ltp — e'^sia 



^ a ;g n (ze lv ) V^e 21 ^ - f 2 j To (35) 
and in the internal region 

k 



(z — eze 2ltp ) Co 



(36) 



Note that the phase factor e 1 ^ should not be absorbed in 
the square root unless its sign is taken into account. 



Fig. 2. A contour - grey-scale plot of the potential of the slice 
which is indicated by the bold white ellipse. 

Typically we deal with concentric isodensities. How- 
ever, we can describe also shifted slices by a proper trans- 
formation to the new origin zo of the shifted slice. We find 
simply 

I = I(z - z ,z - z ) . (37) 

Obviously Eq. (35-36) and Eq. (37) can be combined to 
describe rotated and shifted slices. In Fig. (3) we show the 
deflection angle due to a rotated slice of constant density 
with major and minor axis a = 2, b = 1 respectively. 
The shading increases with increasing norm of deflection. 
The iso- lines show the homoeoidal symmetry (with same 
axial ratio as the slice) inside and the confocal symmetry 
outside the slice. 

Potentials and derivatives of rotated slices can be 
found if the related equations are transformed as shown 
in this section. 

4. Elliptical models for numerical simulations 

The density distributions of gravitational lenses are not di- 
rectly measurable. However, making special assumptions 
about the mass-to-light variation or taking other sources 
of information (velocity dispersion, theoretical models . . . ) 
projected, elliptical density distributions are appropriate 
for modeling elliptical lenses, i. e. the data obtained from 
observation can be used as direct input to the model. Ana- 
lyzing elliptical galaxies typically yields a set of n elliptical 
isophotes (isodensities) a, with elliptical radius r-^i, their 
elliptical parameter e 8 - , their foci /; and the position angle 



Fig. 3. The norm of the deflection angle for a rotated (ip = 30°) 
slice with major and minor axis a = 2, b = 1 is shown. The 
shading increases with the norm (dark ~ high value of the 
norm). The iso-lines show line of constant norm (white) and 
argument (black) respectively. Note the homoeoidal symmetry 
inside the slice and the confocal symmetry outside. The argu- 
ment of the complex deflection is a linear function inside and 
turns to a hyperbola outside. 

of their major axis (twist angle) <fi. We find the deflec- 
tion 7 tot at a location z by a summation over all slices of 
thickness A<7 8 ' = a, — <Ji-\. 

Jtot = Jtot + Jtot (gg) 

All internal contributions are found if 

//Re(ze-^')\ 2 /Im(ze-^')\ 2 , , 

from Eq.(36) 

4 ot =X)/i„(A<7 i ,e i ,p i ) (40) 

i 

and all external contributions for these slices for which 
r(z) > r~Ei from Eq.(35) 

Ilf =X)/ex(A<7 i ,/ i , ¥ » i ) (41) 

i 

Derivatives and potentials can be found respectively. 



Let us take an elliptical density distribution described by 
the set of consistent functions introduced before and pos- 
sibly bounded by an ellipse with radius r max 



e 



c(te) 
e(r E ) 

f = f( r E) 

and in addition 

/ = /M 



(42) 



In the limit of infinitely thin slices Ac — ► da introducing 
the derivative of a, we replace 



d<r(r E ) = <r'(r E )dr E 



(43) 



to find the limit of Eq.(38) using condition Eq.(39) and 
r m ; n = (compare note after Eq. (9)) 



r(z) 



o 2e(r E ) L 

s ^(''E) sign ^g^('-E)^ ^2 e 2i^(r E ) _ j( rE )2 



(44) 



r(z) 



X dr E a'(r E ) (z - e(r E )ze 2l ^ rE )) (45) 

(46) 



and 

/tot = Jtot + pot + / ((T(rmax)) 



where /(c(r max )) is the contribution of the slice of thick- 
ness <r(r max ) which is obviously zero if we deal with un- 
bounded densities which decrease to zero at infinity. This 
term could be omitted if the proper derivatives of the den- 
sity are invoked at its steps but for practical purposes we 
keep it and take simply the proper one-sided derivatives. 
For bounded distributions we have to take the internal or 
external solution (Eq. (17) or Eq. (15))for that single slice 
if the location of z is inside or outside of the the distri- 
bution, respectively. For the latter case we note 7'°' = 
since a' = in that region. The remaining integral I f e ° f ' is 
then to be performed from r = to r max for the same rea- 
son. Note also that an integration from r m ; n = to r max 
changes the direction of integration (with respect to a) 
which introduces an extra minus sign since we summed up 
all slices from bottom to top. This choice is convenient for 
concentric, monotonically decreasing distributions. How- 
ever, our discussion is evidently true for arbitrary distri- 
butions if the boundaries are properly taken into account. 
See Fig. (4) for the geometrical meaning of all contribu- 
tions. 

To find the proper derivatives of Eqs. (44-46) one can 
apply the rules given in section (3.3). 

5.1. Recovering BK 

BK's result for homoeoidal symmetry is easily re- 
evaluated from Eq. (46) if we take e and <p as constants. 



Fig. 4. The different contributions for the deflection of a den- 
sity distribution as given by Eq.(46) at a location z are shown. 
The internal part gives the contribution of all slices in which 
z lies. The external part then gives all contributions to which 
z lies external and the remaining part is due to the thick slice 
given by the cut-off at r max 

For a position inside the distribution (r(z) < r max ) Equa- 
tion (44-45) yield (restricting to the right half plane for 
brevity) 
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Integration by parts gives 
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-{z - ez) [cr(r max ) - cr{r(z))] 



(48) 



(49) 



Adding the contribution of /(c(r max )) we see that all 
terms cancel (using the same argument as in the deriva- 
tion of Eq. (17)), except the remaining part which is just 
Eq. (11) as obtained by Bourassa and Kantowski. Similar 
arguments apply for a location outside the distribution if 
the comments after Eq. (46) are taken into account. 

5.2. Confocal analytical models 

For models with confocal symmetry the foci of each iso- 
density are at the same place by definition, i.e. / and <p are 



With r mm > f/2 we obtain for Eq. (45) and Eq. (44) 
(omitting Bray's correction for brevity) two expressions 
which are related to a density distribution which is con- 
stant from zero to r m ; n and then follows o-(r E ). The vary- 
ing remaining central part could be described by the gen- 
eral method explained above or for special (homoeoidal) 
cases by the BK or KK approach. 
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(51) 



These integrals lead obviously to closed solutions for many 
density distributions. For a pseudo singular isothermal 
sphere 



' psi 



2<r 
2r E -/ 



r E > 



/ 



(52) 



find, e.g. 
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6. Conclusion 

We showed how to derive the lensing properties of an el- 
liptical slice of constant density from previous results of 
BK or Paper I. These slices are taken as basic elements 
to construct numerically or analytically models for the 
deflection due to surface density distributions of general 
elliptical symmetry. The results can be applied to observa- 
tional data, assuming a constant mass-to-light ratio, tak- 
ing the properties of the isodensities (density, ellipticity, 
twist angle) as input. For analytical purposes general one- 
dimensional integrals are derived for the deflection and 
applied (as an example for a closed solution) to the two- 
dimensional pseudo singular isothermal sphere with con- 
focal symmetry. 
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